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We consider three-loop /3-function corrections to the sparticle spectrum in the MSSM, 
with particular emphasis on Snowmass Benchmark points. The three loop running has 
little effect on the weakly interacting particle spectrum, but for the squark masses the 
effect can be comparable to, or greater than, that of two loop running. We extend the 
analysis to the semi-perturbative unification scenario, where the impact of the three loop 
corrections becomes even more dramatic. 
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1. Introduction 



Softly broken supersymmetry remains a well motivated and popular playground for 
Beyond the Standard Model practitioners. Calculations of sparticle spectra resulting from 
given assumptions about the underlying theory have become increasingly refined, with 
several public programs available that incorporate two-loop Renormalisation Group Equa- 
tions (RGEs) and one-loop radiative corrections. For a recent comparison of the output of 
some of these programs see the paper by AUanach, Kraml and Porod (AKP) [|l|. Two loop 
corrections are also available to, for example, the effective potential ^j. One area where 
there has been considerable progress in the last few years is in the calculation of the RGE 
/3-functions. For an arbitrary supersymmetric theory the chiral supermultiplet anomalous 
dimension 7 is known to three loops and the gauge /3-function(s) (3g to four loops |^. 
For the MSSM, both (3g. and the various 7s were given through three loops in Ref. 0. 

The /3-functions in a general theory and the MSSM for the "standard " soft breaking 
terms were given to two loops in Ref. H and for the "non-standard" terms also to two 



loops in Ref. [0. Recently, however, it has been realised^] -||T3] that in the case of the 
"standard" terms it is possible to express the associated /3-functions exactly in terms of 
simple differential operators acting on l3g. and 7. It is therefore a straightforward matter 
to derive the three- loop "standard" soft /3-f unctions for the MSSM and variations thereof. 
In this paper we present three-loop running results for the MSSM with the addition of 
and Tiio sets of SU5 5(5) and 10(10) representations respectively. A motive for grouping 
the additional matter in this way is that complete SU5 representations do not (at one loop) 
change the prediction of sin^ ^VK (or alternatively of g^{Mz)) that follows from imposing 
fi'1,2,3 gauge unification. Also unchanged at one loop is the gauge unification scale, Mx; 
but at higher loops this scale increases and can approach the string scale. (For a recent 
account of unification at the string scale via the addition of incomplete SU^ multiplets, 
see Ref. [0.) 

At three loops each soft /3-function has many terms and some of the coefficients are 
quite large; given this it is worthwhile checking whether perturbation theory remains good 
in the MSSM. It is generally believed that the perturbation series for QFT /3-functions are 
asymptotic in nature. The exact Pg in the NSVZ scheme []T^ for a pure (no matter) = 1 
theory is clearly an exception, but in the presence of matter the perturbation series for 7 
(and hence also the one for (3g) is probably asymptotic (for a discussion see Ref. [|13|)- It 



is interesting that even for 77,5 = nio = we find that, for the squark masses, three loop 
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running corrections are typically larger than the two loop ones. We show explicitly how 
the three loop corrections affect the spectrum for some of the Snowmass benchmark points 

(sps) H. 

Another motive for the calculation and in particular for extending it to n5,nio 7^ 
is that it enables us to explore the phenomenon of "semi-perturbative unification" as de- 
scribed by Kolda and March-Russell (KMR)|T^]. In this scenario, the gauge couplings 
increase at high energies but do not quite reach a Landau pole at gauge unification (con- 
trasting with the non-perturbative unification of Maiani et al[]T^, where the Landau pole 
occurs at Mx)- It is then possible to argue that there is a regime where perturbation 
theory remains reliable, but the resulting physics differs markedly from that obtained in 
the MSSM easel 

Our calculations improve on those of KMR by including one loop threshold corrections 
and the complete three loop running corrections; we check that (for 77,5 = nio = 0) our 
results are consistent with those presented for the SPS points in Ref. [0 and Ref. [0. While 
we provide more precise results, we support the conclusions of KMR by demonstrating that 
^^55 ^10 7^ can lead to changes which are in some cases non-negligible, but are consistent 
with perturbation theory (modulo issues associated with the squark masses which we will 
discuss later), and can be readily distinguished from the MSSM. 



2. The Soft Beta functions 

The procedure for calculating the soft f3 functions from (3g. and 7 is described in 
Ref. 0. The only subtlety relates to the X-function which arises in the soft scalar mass 
/3-function; expressions for the leading and sub-leading contributions to this appear in 
Ref. 0. Armed with these results it is straightforward to calculate the three-loop MSSM 
soft /9-functions from the three- loop expressions for the f3g., 7 given in Ref. 0. We have 
generalised this whole calculation to n5,nio 7^ 0. The resulting expressions are very 
unwieldy; as an example we give the one, two and three loop results for /3^2 , in the 

Qt 

approximation that we retain only (73 and the top quark Yukawa coupling At (in what 
follows we denote the third generation squarks as Qt,t^,b'^, the first or second generation 
squarks as Qu,u'^-,d^-, and we suppress IGtt^ loop factors. ): 



For other related work on unification at strong coupling see Ref. | 17 
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=2A2(E, + A?)-8(^^fMf + |^7iM| + |^|M|) (2.1a) 

Qt 

= -20A4(E, + 2^2) + WglMiin^ + 3nio - f ) 

Qt 

+ fali^ml^ + ml + ml + (mo + 2){ml. + 2m|J + (m + 2)ml) {2.1b) 
/32 = [(1280fc + ^ + 16ni + (^ + f^fc)(n5 + 3mo) 

Qt 

+ 96nion5 + 144n?o)M| + (^ - ^(ns + 3nio))(m?c + m^ + 2m|J 
+ (2mJ^ + ml){^ - f ns + f nio - fn^nio - IQnjo) 
+ mU^ + - 32nio - IQn^nio - fnjM 

- [(288 +^k + 48(n5 + 3nio))M| - (192 + if^A; + 32(n5 + 3nio))AtMs 
+ (272^ + 176 + 8(n5 + 3nio))(E, + A',)]X',gl 

+ (if + 32A;) [Mi - 2AtMs + + 2^^)] X^gj + {6k + 90) (Et + 3A^^)Xl (2.1c) 

where /c = 6^(3), and Et = mg^ +m| + m|c . For this special case, and also with = nio = 
0, the three loop result, Eq. ( |2.1c|) , was given in Ref. fl^, except that in the corresponding 
expressions in this reference the squark masses of different generations are not clearly 
distinguished (as they must be since the third generation evolves differently from the other 
two). Complete results for the three loop /3-functions including all three gauge couplings 
and Ug X Ug Yukawa matrices may be obtained by application to the authors. 

Note that in our analysis we do not include "tadpole" contributions, corresponding 
to renormalisation of the Fayet-Iliopoulos (FI) L>-term. These contributions are not ex- 
pressible exactly in terms of (3g-, 7; for a discussion, and three loop results for the MSSM, 
see Ref. . For universal boundary conditions, the FI term is very small at low energies 
if it is zero at gauge unification; in this paper we restrict ourselves to universal boundary 
conditions and ignore these contributions. 

3. The Running Analysis 

In this section we examine the effect of the three loop corrections on the standard run- 
ning analysis. We will focus on the standard treatment with universal boundary conditions 
at gauge unification, often termed CMSSM or MSUGRA. Thus we assume that at Mx 

we have universal soft scalar masses (mo), gaugino masses (mi) and A-parameters {A), 

2 

and work in the third-generation-only Yukawa coupling approximation. This is for ease 
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of comparison with existing results rather than because we find the scenario particularly 
compelling. In the MSSM the corrections to the dimensionless coupling running analysis 
due to two and three loop corrections are comparatively small[^. As emphasised by KMR, 



however, this becomes less true for 77,5, nio 7^ 0. In particular, for 77,5 + 3n 
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6, /3, 







at one loop, so that we need to consider at least two loop corrections, and also three loops 
to verify that we remain within the perturbative domain. 

We will present results as a function of nio (for 715 = 0); the dependence of physical 
quantities on (for nio = 0) is qualitatively similar though not identical. As remarked 
by KMR, the mass scale of these additional multiplets being unknown it makes sense to 
parametrise their effects by taking n^^nio to be continuous variables. 



.a^.a^ vs 1: for SPS1 a 
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Fig.l: Gauge coupling unification for nio = 1.7. Solid, dashed, and 
dotted lines correspond to ai,a2,as respectively. 

In Fig. 1 we show the evolution of the gauge couplings = gf/{47r) for nio = 1-7, 
using three loop /9-functions for all couplings. The couplings are plotted against r = 
^ln{Q/Mz)', evidently we are still in the perturbative regime. The input parameters at 



Mz correspond to a typical supersymmetric mass spectrum; specifically, the Benchmark 
point SPSla, the details of which will be given later. We adjust input parameters according 
to the supersymmetric spectrum in order to account for threshold corrections in the manner 
of Ref. [P i, and run up from Mz using the full supersymmetric /3-functions; thus the 



input values for the gauge couplings depend on the sparticle spectrum, and are determined 
iteratively. We have set nio = 1.7 because around this value the electroweak vacuum 
fails for the SPSla input parameters; at this value it is interesting that (see Fig. 1) the 
(small) one-loop contribution to /3q,3 is almost precisely cancelled by the two and three- 
loop ones. As already mentioned, if we have gauge coupling unification then this (and 
the scale at which it occurs) is unaffected by taking n^^nio 7^ in the one loop running 
approximation. This ceases to be true for two and three loop running, and we must 
assume the existence of quite large GUT-scale threshold corrections to ensure unification. 
The unification scale corresponding to Fig. 1 (defined as where ai and a2 meet) is Mu ~ 
1 X lO^^GeV, significantly higher than in the MSSM. 

Turning to the soft parameters, note that there are some large coefficients in the ex- 
pression for P^^l in Eq. ( |2.1c] ); the coefficient of the M^g^ term, for example is O(IO^), 

Qt 

even in the MSSM when ns = tiiq = 0. For weakly-interacting particles and the gluinos, 
the three loop effects are quite small for zero or small 77,5, nio; but for the squark masses 
the three loop /3-function coefficients are (at Mz) typically (while smaller than the corre- 
sponding 1-loop coefficients) larger than the corresponding two loop coefficients (even at 

^5,10 = 0) if mi > mo, as will be so, in fact, in the cases we shall present. One might well 
2 

be tempted to interpret this as evidence for the asymptotic nature of the /5-function series, 
as we mentioned in the Introduction. We will see the effects of this in the next section. 
As mentioned above, we adjust the input dimensionless parameters to accommodate 



threshold corrections [21]. We also incorporate the one-loop radiative corrections from this 



reference. For the input top mass at the weak scale we use [p2| : 



mt(Q) =mt , \l - —ih - ?,L) - a\( 0.538 - ^ + 
+ 0(0^3) -t- electroweak, sparticle contributions 



43L 3L2- 

' (3.1) 



^ In the first line of Eq. 37 of Ref. tlie first term in tlie square bracket sliould read 

— (m^^ + m^^)Bo{m^^,m^_^,0): i.e. it should have a minus sign. The corresponding exact result in 
Eq. D49 is correct, however. 
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where L = ln[mt(Q)^/(5^]. This formula is identical to one from the up-to-date version of 
AUanach et al|I[]. 

For definiteness we analyse in detail particular SPS benchmark points[n|, and verify 
that our results (for = nio = 0) are in accord with those obtained with existing 
computational tools [Q, |1T8[]. 



3.1. Benchmark point SPS la 

This point is a "typical" point in MSUGRA parameter space, with mo = lOOGeV, 
mi = 250GeV, Aq = — lOOGeV, tan/3 = 10 and /i > 0. In Table 1 we compare our results 

2 

for a selection of sparticle masses (at = nio = 0) with the spread of results quoted in 
AKP (note our convention that the predominantly R-handed top squark is ^2)- 



mass 


Hoop 


2loops 


"iloops 


AKP 


9 


630 


615 


612 


594 - 626 


i2 


404 


403 


395 


379 - 410 


ul 


571 


563 


555 


536 - 570 


ur 


551 


547 


538 


520 - 569 


LSP 


105 


97 


97 


96.4 - 97.6 



Table 1: Sparticle masses (in GeV) for the SPSla point 



We would expect our two loop results to correspond most closely to AKP and we 
see that they are indeed consistent. The effect of inclusion of three loop running is never 
greater than 2%; note, however, that the shift caused by three loop running effects is 
comparable for ul and larger for t2,UR than that produced by two loop running effects. 
In Fig. 2 we plot the ratio of the light stop mass to the gluino mass as a function of tiiq 
(for ns = 0), and we see that both two and three loop effects increase dramatically as tiio 
increases. (NB the increase in the squark/gluino mass ratio with n^^iQ observed by KMR 
applies to the squarks of the first two generations). 

In Fig. 3 we plot the ratio of the LSP mass to the gluino mass as a function of nio (for 
77-5 = 0). Here we see that the impact of the three loop running corrections is less marked 
but still appreciable at large niQ. 

7 



0.75 



Light stop/gluino mass ratio vs n for SPS1 a 
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Fig. 2: Plot of the light stop/gluino mass ratio against Uiq for SPSla. 
Solid, dashed and dotted lines correspond to one, two and three loop 
running respectively. 
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Fig. 3: Plot of the LSP/gluino mass ratio against nio for SPSla. Solid, 
dashed and dotted lines correspond to one, two and three loop running 
respectively. 
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3. 2. Benchmark point SPS 5 



This point differs from the previous one in having a large value of the ^-parameter: 
with mo = 150GeV, mi = 300GeV, Aq — — ITeV, tan/3 = 5 and > 0. The contributions 

2 

of HiAq to the off-diagonal term in the stop mass matrix have the same sign, and the 
magnitude of is large, resulting in a light stop. For this point we obtain in the MSSM 
(with 77,5 = 77-10 = 0) the results shown in Table 2. 



mass 


Hoop 


2loops 


Sloops 


AKP 


9 


745 


731 


729 


705 - 730 




236 


250 


231 


232 - 248 


ul 


682 


674 


666 


642 - 681 


ur 


657 


655 


645 


622 - 681 


LSP 


128 


120 


120 


118.7-121.1 



Table 2: Sparticle masses (in GeV) for the SPS5 point 



Once again the inclusion of three loop effects causes mass shifts of less than 2%, except 
for the light stop where there is an effect of about 8%. The light stop mass comes from a 
2x2 matrix with large off-diagonal entries; the large three loop shift is caused essentially 
by the changes in all the entries due to the comparatively large contributions to the three 
loop soft /3-functions, as we described earlier. The mass of the light stop is also very 
sensitive indeed to the input mt{Mz), which in turn depends on the sparticle spectrum 
and the input top pole mass (which we have taken to be 174.3GeV). However (for a given 
™tpoie)' ™t(-^z) changes very little when we include the 3 loop corrections. Note that 
our two loop result is slightly above the range obtained by AKP. It is also very sensitive 
to 77-5,77-10; in Figure 4 we again plot the light stop/gluino mass ratio against uiq: this 
time the electroweak vacuum fails for 77io ~ 0.3. Remarkably, the three loop and two loop 
corrections cancel almost exactly near riio — 0, while at uiq ~ 0.3 the two loop corrections 
are rather small. 
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Light stop/gluino mass ratio vs n^^ for SPSS 
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Fig. 4: Plot of the light stop/gluino mass ratio against nio for SPSS. 
Solid, dashed and dotted lines correspond to one, two and three loop 
running respectively. 



LSP/giuino mass ratio vs n.^^ for SPS5 
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Fig. 5: Plot of the LSP/gluino mass ratio against nio for SPSS. Solid, 
dashed and dotted lines correspond to one, two and three loop running 
respectively. 
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In Fig. 5 we plot the ratio of the LSP mass to the gluino mass for SPSS as a function 
of nio (for 77-5 = 0). Here we see that the impact of the three loop running corrections is 
less marked but still appreciable as nio increases. 

4. Conclusions 

We have extended typical detailed running coupling analyses for a selection of MSSM 
SPS benchmark points to incorporate three loop /9-function corrections for the running 
masses and couplings. Generally speaking the effect of the three loop running corrections 
is at most 2% and of the same size or smaller than that of the two loop corrections, 
except for squark masses where it can be larger; simply because the three-loop /3-function 
coefRcients are larger than the two-loop ones. For the light stop mass for the SPS5 point, 
we see an 8% effect; this happens because this mass results from the diagonalisation of a 
matrix with large off diagonal entries which all change as described above. We have also 
performed the same analysis for the MSSM extended to incorporate additional matter in 
the form of SU^ 5 and 10 representations. As the amount of such matter is increased the 
effect of two and three loop corrections becomes more dramatic, as the one-loop /^-function 
for Q!3 decreases. 

Given some detail of the sparticle spectrum, it will be comparatively easy to dis- 
tinguish, for example, the CMSSM and AMSB scenarios; however in the context of the 
former, disentangling the possible impact of additional matter and the effect of radiative 
corrections will be more difficult. 

One may expect in general that two loop threshold/pole mass corrections will be 
competitive with the three loop running corrections that we have described, and so for 
accurate predictions one should include both. At this level one is also sensitive to the 
experimental uncertainty in mt (for the light stop sometimes very sensitive, as we have 
described) and the strong coupling as{Mz)- It is feasible that by the time sparticles are 
discovered complete two loop threshold corrections will be available, and that both these 
uncertainties will also be reduced, so that significantly more accurate sparticle spectrum 
predictions will be possible. It appears, however, from the apparently asymptotic nature 
of the squark mass /^-functions that squark mass predictions with an accuracy greater than 
around 2% will not be possible using perturbation theory. 
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